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Efficient laser ion acceleration requires high laser intensities, which can only be obtained by tightly 
focusing laser radiation. In the radiation pressure acceleration regime, where the tightly focused 
laser driver leads to the appearance of the fundamental limit for the max;imum attainable ion energy, 
this limit corresponds to the laser pulse group velocity as well as to another limit connected with 
the transverse expansion of the accelerated foil and consequent onset of the foil transparency. These 
limits can be relaxed by using composite targets, consisting of a thin foil followed by a near critical 
density slab. Such targets provide guiding of a laser pulse inside a self-generated channel and 
background electrons, being snowplowed by the pulse, compensate for the transverse expansion. 
The use of composite targets results in a significant increase in maximum ion energy, compared to 
a single foil target case. 
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I. INTRODUCTION 

The laser-driven acceleration of charged particles is the main focus of many existing, being constructed, and projected 
high-intensity laser facilities |l|43|. In particular, the acceleration of ions has attracted a significant attention both 
theoretically and experimentally W, due to the many potential applications of laser-accelerated ion beams, including 
but not limited to, fast ignition Q, hadron therapy radiography of dense targets Q, injection into conventional 
accelerators §, and nuclear physics studies [§|. 

Several years ago Particle-in-Cell (PIC) based computer simulations predicted that several hundred MeV or even 
GeV ion beams can be produced in laser-matter interactions using modern multi-TW or PW laser systems M- 
Recently, several new acceleration schemes were proposed that offer the possibility of generation multi-GeV or TeV 
proton beams EH- However such schemes would require next generation lasers with even higher peak intensity. 
There are several basic laser ion acceleration mechanisms that has been discussed in the literature: (i) Target Normal 
Sheath Acceleration (TNSA) [l^, (ii) Goulomb Explosion (GE) [l3|, (iii) Radiation Pressure Acceleration (RPA) [l5| . 
(iv) Magnetic Vortex Acceleration (MVA) [T6l - [T8l |. and (v) the Shock Wave Acceleration (SWA) [II- 

TNSA is usually realized in the case of relatively thick targets and poor laser contrast, producing ion beams with 
Maxwellian-type energy spectra. In TNSA the laser launches electrons from the front of the target all the way through 
the target. As they leave the target from the back, those electrons establish a quasistatic longitudinal electric field 
that accelerates ions from the rear side of the target. Most of the experimental results on proton acceleration can be 
attributed to the TNSA scheme, with the peak proton energy around 70 MeV [^ . 

The GE acceleration works for high intensity pulses that are able to remove all the electrons from the irradiated 
spot, ideally without disturbing the ion core. The remaining ion core experiences GE due to the noncompensated 
positive charge. Usually a two layer (high Z/low Z) target is used to enhance the properties of GE accelerated ions. 
The low-Z ions are accelerated in the field of expanding high-Z layer producing a quasimonoenergetic spectrum, which 
is required for almost all laser driven ion beam applications. However extremely high intensity and high laser contrast, 
not only at the nanosecond scale but also at the picosecond scale make this scheme very challenging for experimental 
realization. 

The RPA mechanism comes into play when the laser is able to push the foil as a whole by its radiation pressure. 
It is the realization of the relativistic receding mirror concept (^ . The role of a mirror is played by an ultra-thin 
solid density foil or by plasma density modulations emerging when the laser interacts with an extended under-critical 
density target, the so-called hole-boring RPA [22|. The laser pulse is reflected by a co-moving mirror. The problem 
of a plane EM wave reflection by a mirror moving with a relativistic velocity was considered by A. Einstein as an 
illustration of the Theory of Special Realtivity [^. The frequency of the reflected radiation is shifted down by a 
factor of 47 ^, where 7 is the Lorentz factor of the mirror. Thus the energy transferred to the mirror is (1 — 
where £l the energy of the laser pulse. For 7^1 almost all laser energy is transferred to the foil, which makes this 
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scheme very attractive in the ultrarelativistic limit [l^, [2^ . Recently, there were several papers published that 

claim the experimental observation of the onset of this regime of laser ion acceleration [2^. 

The MVA regime is different from the other three regimes, and is realized when the laser interacts with a Near 
Critical Density (NCD) target (TNSA, CE, and RPA rely on the interaction with solid density ultra-thin foils). In 
MVA the laser pulse propagates in NCD plasma, it makes a channel in both electron and ion density. When the laser 
pulse exits the plasma it establishes strong longitudinal electric fields that can accelerate the ions which are pinched 
by the electrons in a thin filament along the laser propagation axis. In the case of long pulses, the acceleration of 
helium atoms up to 40 MeV from underdense plasma was observed at the VULCAN laser and the acceleration 
of protons up to 50 MeV at Omega EP [2^. Also the experiments with short pulses and cluster jets show that 10-20 
MeV per nucleon ions can be generated via such interaction [2^ . Like MVA the SWA regime employes an underdense 
plasma, where the laser launches a shock wave inside the target. The ions in the bulk of the target are reflected 
by this shock wave at twice the velocity of the shock. Ideally such mechanism should produce monoenergetic ions, 
since the reflection happens after the shock is detached from the laser pulse. However, the accelerated ion beam will 
interact with the background plasma and will be subject to two-stream and filamentation instabilities the latter 
determining the width of the accelerated ion spectra. 

There are several mechanisms that, through either modification or combination of some of the basic regimes, 
enhance the maximum ion energy, number of accelerated ions, or improve their spectrum. For example, the Burn-out- 
Afterburner (BOA) [^, which employes an enhanced TNSA, and Directed Coulomb Explosion (DCE) [s^, which 
is the combination of RPA and CE, are such composite mechanisms. There was also an extensive study connected 
with ion acceleration during the induced relativistic transparency of the target (s^. Also the use of composite targets 
(low density/high density) was proposed in a number of papers to either inject the ions into accelerating fields or to 
enhance the interaction of the laser pulse with the high density part of the target [Hill HE [Hi. 

In this paper we propose to enhance the maximum attainable energy of ions accelerated via RPA by utilizing 
composite targets, consisting of a solid density ultra-thin foil placed in front of the NCD slab. The RPA accelerated 
ions are subject to a number of factors limiting their maximum energy. Typically RPA requires high laser intensities, 
which can only be obtained by tightly focusing laser radiation. In this case the finite spot size effects are important 
(H |. Tightly focused laser pulses have a group velocity less than the speed of light in vacuum, which would limit 
the maximum attainable ion energy to the energy corresponding to the group velocity [H . [H | . In the framework of 
the receding relativistic mirror concept, the energy transferred from the pulse propagating with the group velocity, 
Pg, smaller than the speed of light in vacuum to the mirror is A£ ~ 27 ^/ 3 (/ 3 g — P)£l- Thus, if /? = /3g, then there is 
no interaction of the laser light with the target, and consequently no energy transfer. Moreover tightly focused laser 
pulses force the foil to expand in the transverse direction during the interaction, making it transparent for radiation 
and effectively terminating the acceleration. In the case of a composite target, the laser, after beginning to accelerate 
the foil, enters the NCD slab and generates a channel in electron and ion density. Inside this self-generated channel 
the laser pulse propagates at subluminal group velocity, Pg < 1, which is a function of the laser pulse power, over 
the depletion length, which is also the acceleration length for the foil. For a tightly focused pulse, the transverse 
expansion of the foil is the main limiting factor for the acceleration. The acceleration distance for the expanding foil 
is substantially smaller than the depletion length in the NCD plasma. Thus the interaction of a laser pulse with a 
composite target leads to the enhancement of the maximum ion energy. 

The paper is organized as follows. In section 2 we review the equations of the foil motion under the action of the 
laser pulse radiation pressure. In section 3 we consider the effect of the transverse expansion of the target. We present 
the results of PIC simulations in section 4 and conclude in section 5. We discuss the properties of the laser pulse 
propagating in the NCD plasma and the energy of ions accelerated via the MVA mechanism in Appendices A and B 
respectively. 


II. RADIATION PRESSURE ACCELERATION 


In this section we review the basic properties of the RPA mechanism. The motion of the foil is modeled by an 
interaction of an EM wave with a mirror described by the equations [HIH,!!!. The laser pulse radiation pressure 
is taken into account as a force proportional to the EM wave momentum, which in turn is proportional to the EM 
wave Pointing vector, S = E/, x Bi/47r, where El and Bl are the electric and magnetic fields of the laser pulse. We 
assume c = 1 throughout the paper. If we consider the RPA by a plane wave of a non-expanding opaque foil in the 
ID approximation |l5l |. then the momentum of the foil is 
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is the laser pulse fluence (incident laser energy per unit area) and ngl is the areal density of the foil, with Ue the 
electron foil density and I the foil thickness. In the ultrarelativistic case the energy of the foil asymptotes to [a 
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The subscript Pg = \ indicates that the laser group velocity is equal to the speed of light in vacuum in this case. 

In a more general case of an expanding target and a laser pulse with group velocity smaller than the speed of light in 
vacuum the equation of motion of the surface element of a mirror in the laboratory frame of reference can be written 
in the following form [s^, [13 : 


dt riel ’ 


( 4 ) 


where the force is proportional to the laser pressure and directed along the normal to the surface element, denoted by 
unit vector v. The light pressure is determined from balancing the fluxes of the incident, reflected, and transmitted 
EM waves in the reference frame, where the particular surface element is at rest. In what follows we will refer to this 
reference frame as the moving frame or M-frame, and to the laboratory frame as the L-frame. In the M-frame the 
force acting on the foil is 

F = (l + |p|2-|rn(S.r.)r.. (5) 


The reflection coefficient, p, the transmission coefficient, r, and the absorption coefficient, a, are related to each other 
through energy conservation: |pp + |rp + jap = 1, and are determined from the solution of the wave equation for 
the EM wave interacting with a thin foil in the M-frame [13 ■ 

The maximum attainable ion energy for RPA by laser pulses having their peak intensity and consequently peak 
fluence on-axis is also achieved by an on-axis surface element of the foil. For this element the laser propagtion 
direction coincides with the vector normal to the surface. For a circularly polarized EM wave with amplitude of the 
vector-potential, Aq, frequency w, and wave vector k, propagating along the x axis we have 

S = cjkAle^, (6) 

where Aq is the amplitude of the EM wave vector potential. Following the results of Ref. [sO], we assume that the 
EM wave travels with a group velocity jdg = k/uj^ which is the case for focused pulses and waves traveling inside some 
external guiding structure. In the boosted frame of reference, moving with velocity /3, we obtain: 


ujk = 




I-/32 


( 7 ) 


Thus the equation of the on-axis surface element motion is 

where time is measured in units of uj~^, and 
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Here a = eAjme is the dimensionless vector-potential, eg = TT{nel/ncrE) is the parameter governing the transparency 
of the thin solid density target (^ . Tier = meW^/47re^ is the critical plasma density, e and me are the electron charge 
and mass respectively, i!) = {t — x/jdg) is the phase, and mi is the ion mass. The reflection coefficient is determined 
from the solution of the wave equation in the rest frame of the foil : 
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^ ^ - eo - 1)^ + 4ag] + («o “ “ 1)} ^ > (H) 

where 7 = (1 and eg = 760 , which accounts for the fact that the parameter eg is not a relativistic invariant 

( 3 ^ . We notice that the condition ap = ep is considered to mark the optimal regime of ion acceleration 0- However 
this is true only in the case of non relativistic ion energies. In the relativistic case, the condition transforms into 
Op = 7 ep [ 13 . In what follows we assume total reflection and set |p| = 1. For /3 = /3g, the right hand side of Eq. (8), 
i.e., the radiation pressure is equal to zero. This will result in the appearance of the maximum attainable velocity 
equal to (3g for the ions accelerated by the RPA of an EM wave with Pg < 1. The equation (jS]) can be solved in 
quadratures. For the sake of brevity we assume that initially the foil was at rest, /3(0) = 0. This yields the following 
expression: 


In 


( 1 -/ 3 ^, + ( 1 -/ 32 ) 1 / 2 ( 1 -/ 32 ) 1 / 2 ^^ 

(/3,-/3)(l + (l-/32)i/2) 
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= Pg{l-Plf/^Kp{t)., 
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t 

where Kp(t) = f ndt'. In order to get an idea of the foil’s behavior as its velocity approaches Pg, we assume that the 
0 

EM field is constant, then Kp(t) = nt and for t —>■ 00 only the term with ln(, 5 g — P) survives. Thus we obtain the 
following asymptotic expression for the ion velocity: 

P = Pg-e-xp[-Pg{\- Plf/'^Kt^. (13) 

From Eq. m one can see that the foil velocity approaches the laser group velocity exponentially, but never exceeds 
it. Thus the ion energy is limited by 7 g = (1 — Pg)~^^'^ and the time needed to accelerate the foil to the energy 
approaching 7^ scales logarithmically with 7: 
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We should mention here that, for finite duration laser pulses, the group velocity limit will manifest itself only for 
l^g=i > lgi i- 6 -i in the case where the maximum accelerated ion energy by the same pulse (but with Pg = 1 ), would 
exceed the energy determined by the group velocity. 

In Fig. 1 we present the numerical solution of Eq. m for a Gaussian pulse with a duration of 27 fs (10 cycles), 
A = 800 nm wavelength, focal spot of wq = 1.35A, which corresponds to an f-number of f /D = 1.5, interacting with 
a 0.25A thick hydrogen foil with an electron density of 40077.^^. The evolution of the maximum ion velocity (Fig. la) 
and maximum ion energy (Fig. lb) are shown for three different values of the averaged laser power (1 PW, 1.8 PW, 
and 3.6 PW) and two values of group velocity Pg = 1 and Pg = 0.986 {jg = 6 , wq = 1.35A, f/D=1.5 0). For P=1 
PW, the evolution of ion velocity (and energy) is very similar for Pg = 1 and Pg < 1 cases, since < 7 g and 

the effects group velocity are small. For P=1.8 PW, one can see a small difference between these two cases, which 
is mainly due to different dephasing between the pulse and the foil lengths. We will address this effect below. For 
P=3.6 PW, > 7 g, which leads to significant differences between Pg = 1 and /3g < 1 cases, as expected from the 

analytical analysis. The Pg < 1 cases are limited by Pg for velocity and 7 g — 1 for energy, while the Pg = 1 cases are 
only limited by the total energy stored in the laser pulse. 

We notice that for P=1.8 PW, the Pg < 1 case results in higher maximum ion velocity (energy). This is due to the 
fact that, when the foil reaches relativistic velocities, the dephasing between the foil and the laser plays an important 
role. The pulse with Pg < 1 will interact with the foil longer, thus accelerating it to higher energies. If we rewrite Eq. 
dS]) in terms of the phase tp = {t — x/Pg), then 


^=nPl{l-P^)^/\l-PPg). 

This equation can also be solved in quadratures: 
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Ip 

where Kp{ip) = f ndtp'. From Eq. (fTHl) we obtain for the surface element velocity: 
0 


fin sin^ R — (1 — fif) cos R 

P = - —^— m — 

1 — Pg cos"^ R 

where 

R = Ki 3 {'il;)Pg{l - + arccos/3g. (18) 

Since 0 < fi < Pg, the parameter R is defined on the interval arccos/3g < R < tt/2. In order to compare the cases 
Pg = 1 and Pg < 1 we assume that Kp{p)) = Kif and write down the expression for the final velocity, in the case 
/3g = 1 , in the following form 


Kip (2 + Kip) 

2 + 2kiP + {kiP)^ 


(19) 


In the case of /3g < 1 we use Eq. m, assume that /3g —>■ 1 , expand the lefthand side around Pg = 1, and solve with 
respect to P, yielding 
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Thus one can see that, due to slower dephasing, the laser pulse with Pg < 1 accelerates the foil to higher energies 
than the laser pulse with Pg = 1. However this consideration is valid only for P < Pg. This dephasing is illustrated 
in Fig. Ic, where the evolution of the value of the field amplitude, ao(x), at the foil is shown. For the particular case 
of P=1.8 PW wee see that Pg < 1 curve has higher values than the Pg = 1 curve for an extended period of time. In 
the case of a 3.6 PW laser, the pulse becomes phase locked with the foil, i.e., it co-propagates with it, but can not 
transfer any amount of energy to the foil. 

From Eq. we can estimate the duration of the pulse for a given laser power for the pulse to become phase 
locked with the foil. We write down an equation for the phase corresponding to P = Pg: 


Kpi^max) - _ ^2)1/2 (2 arccos/3g^ , 

which for Kp{ip) = Kip, i.e., for constant EM field, is reduced to 

^max = ^^2(1 _ ^2)1/2 (3 " ^rCCOS/Sg) , 


( 21 ) 


( 22 ) 


or for Pg —^ 1, ipmax ~ (^/2)(7g/K). The Eq. (21) indicates the minimal pulse duration for a given intensity that is 
needed to accelerate a foil up to the velocity equal to the group velocity of the laser pulse, which basically states a 
well known fact that, in the ultra relativistic case, almost all the laser energy is transferred to the ions [l^. Moreover 
Eq. (jir)|) is defined on the interval 0 < ip < ipmax- Since we can rewrite k as k = the maximum phase is 


Ipmax — „ mix (2^) 

If (7r/2)7g > then ipmax > t and the foil can not be accelerated to 7 g, since there is not enough energy in the 

pulse. In this case the difference in maximum energy evolution in the cases of /3g = 1 and /3g < 1 should be small. 
For Ipmax < T the difference should be rather large, since for /3g < 1 the laser pulse would be able to accelerate ions 
only up to 7 g — 1 , while in the case of ,5g = 1 the ion maximum energy is determined by the total energy of the 
pulse. The dependence of the maximum ion energy evolution on the duration of the pulse for two cases of /3g = 1 and 
/3g < 1 is illustrated in Fig. Id. Here the laser pulse is chosen to have a rectangular longitudinal profile. The curves 
corresponding to 15 fs laser pulse duration are indistinguishable for /3g = 1 and Pg < 1, which indicates that the group 
velocity effects are almost negligible. In the 30 fs case there is a slight difference, which is due to the dephasing of the 
laser pulse and the foil, and 60 fs and 120 fs cases demonstrate a significant difference. The /3g < 1 curves are limited 
by ( 7 g — 1), while Pg = 1 curves go beyond this value. We can also see that Pg < 1 curves for 60 fs and 120 fs are 
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FIG. 1: (a) The dependences of the foil velocity on time for three values of the averaged laser pulse power: (1) 1 PW, (2) 1.8 
PW, and (3) 3.6 PW for the cases of j3g — 0.986 (solid curves) and Pg = 1 (dashed curves), with the laser pulse duration 30 
fs. The horizontal line marks pg. (b) The dependence of the foil energy on time for three values of the laser pulse power: (1) 
1 PW, (2) 1.8 PW, and (3) 3.6 PW for the cases of Pg = 0.986 (solid curves) and Pg = 1 (dashed curves), with the laser pulse 
duration 30 fs. The horizontal line marks 7 g — 1. (c) The evolution of the laser pulse EM field amplitude at the foil for three 
values of laser power: (1) 1 PW, (2) 1.8 PW, and (3) 3.6 PW for the cases of Pg — 0.986 (solid curves) and pg = 1 (dashed 
curves), (d) The dependence of the foil energy on time for four values of the laser pulse duration: (1) 15 fs, (2) 30 fs, (3) 60 
fs, and (4) 120 fs for the cases of /3g < 1 (solid curves) and /3g = 1 (dashed curves), with the laser pulse power 1.8 TW. The 
horizontal line marks 7 g — 1. The asterisk mark the {Pg = l)-curves. The density of the foil is We = AOOricr and the thickness 
is I = 0.25A. 


almost indistinguishable, since only the At = ipmax < t part of the pulse is able to accelerate the ions. There is no 
such similarity between the 60 fs and 120 fs cases, when /3g = 1 . 

The limiting factor oi Pg <1 can also be described in terms of the phase evolution. It is governed by the following 
equation 


which can be solved in quadratures: 


^ ^ Pg - PW 

dt Pg 
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If t —>■ oo, then depending on the ratio of and 7 g we have the following estimates for the phase evolution, (i) 

When 7 ™“]^ < 7 g, the acceleration processes for Pg = \ and Pg <1 are indistinguishable, R <C 7 r/ 2 : 
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(26) 


(ii) In the opposite case of 7 ™“]^ > 7 g, the acceleration process is greatly affected by the existence of the fundamental 
limit on the maximum attainable ion energy, resulting in a different equation for the phase evolution: 
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4max 




• exp 



(27) 


Here we see that ip approaches 4max but never exceeds it. Thus a laser pulse with duration r > 4max will accelerate 
the foil up to /3 = /3g, but the part of the pulse characterized by the phase values 4max < 4 t will not be able to 
reach the foil and will not take a part in the acceleration process. 
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FIG. 2: The dependence of the phase tp on time for a 1 PW, f/D=1.5 laser pulse for four durations of the pulse: 13 fs (1), 27 
fs (2), 54 fs (3), and 108 fs (4). The maximum phase from Eq. (I23II is plotted as a dashed horizontal line. The density of the 
foil is Tie = 400ncr and the thickness is I = 0.25A. 


In Fig. 2 we demonstrate the dependence of the phase on time. For < jg (curves I and 2) the phase tends 

to infinity as the time tends to infinity, which means that the pulse has accelerated the foil to some energy and is 
no longer interacting with the foil. In the case 7 ™“^^ > 7 g (curves 3 and 4) the phase tends to 'ipmax, he., the pulse 
becomes phase locked with the foil. When phase locked, accceleration has stopped, since the laser is no longer able 
to transfer energy to the foil. Thus the laser pulse and the foil co-propagate without interaction and ip = ipmax 
corresponds to the foil reaching the velocity equal to the group velocity of the pulse. 

Thus, following the approach of Ref. [l2| we showed that there exists a fundamental limit on the maximum 
attainable ion energy in the RPA regime and it is equal to the 7 -factor of the driving laser pulse. This limit greatly 
affects the acceleration process for the pulses with 7 ^“fi > 7 g. We also identified a minimum laser duration needed 
to reach the maximum accelerated ion energy for fixed laser intensity, II'. 


_ TT-fgUel 

'^min — ^ j 


(28) 


III. TRANSVERSE EXPANSION OF THE TARGET 


A group velocity smaller than the vacuum light speed appears naturally in the case of tightly focused laser pulses 
[I|, which diverge rather quickly after passing through the focus. Assuming that this divergence forces the irradiated 
part of the foil to expand, following the increase of the laser spot size, we study how the transverse expansion 
of the target limits the maximum attainable ion energy during the RPA, and whether this limitation dominates 
over the fundamental effects of the group velocity. Assuming that the field of the pulse can be given by the paraxial 
approximation, characterized by the laser pulse waist at focus, wq and the Rayleigh length = tticq/A, the evolution 

of the laser pulse waist as it travels away from focus is w{x) = wq [l -I- {x/Ln)^] , the amplitude of the field scales 

with the distance from the focus as a{x) = ao [l -I- {x/Lji)'^'\ and the group velocity is /3g ~ 1 — 1/k^WQ [i3|. Since 
we are interested in the maximum ion energy, we consider RPA of an on-axis element of the foil. The intensity profile 
near the axis can be approximated by an expanding spherical cup with curvature radius equal to the laser waist, w{x). 
The on-axis element of the foil can also be approximated by an expanding spherical cup with the curvature w{x) and 
areal density equal to rie/ = [l -|- [x/LrY] and ee{x) = ee( 0 ) [l -|- {x/Lr)Y ■ Substituting the field and areal 
density into Eq. ([5]), we see that the right hand side of Eq. ([5]) depends on the distance from the focus only through 
the reflection coefficient [^ : 


(Xl 
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(31) 
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FIG. 3: (a) The evolution of the maximum ion energy for two values of f-number: f/D=l (blue curves, 1, 2, and 3) and f/D=1.5 
(magenta curves, 1*, 2*, and 3*). The solid curves (1 and 1*) are the solutions of Eq. (I29II with laser divergence and target 
expansion taken into account. The dotted curves (2 and 2*) are the solutions of Eq. ([51) with jSg = 1 and dashed curves (3 and 
3*) are the solutions of Eq. m with Pg < 1. (b) The evolution of EM field amplitude at the foil corresponding to the curves 
in Fig. 3a. (c) The evolution of the maximum ion energy for f/D=l (blue curves, 1 and 1*), f/D=1.5 (magenta curves, 2 and 
2*), and f/D=2 (red curves, 3 and 3*). The curves 1, 2, and 3 are the solutions of Eq. (1291) . while the curves 1*, 2*, and 3* 
correspond to the guided case, where the laser is not diffracting and the reflection coefficient is equal to one. (d) The evolution 
of EM field amplitude at the foil corresponding to the curves in Fig. 3c. Laser pulse power is 1.8 PW, duration is 30 fs, the 
foil thickness is 0.25A, and density is Tie = dOOricr. 


In what follows we solve Eq. numerically, taking into account transverse expansion of the foil and laser pulse 
divergence. The comparison between the cases with the transverse expansion taken into account and without is shown 
in Fig. 3 for a 1.8 PW laser pulse interacting with a 0.25A thick hydrogen foil with the density of Ug = 40077.^ for 
three values of the f-number: f /D = 1, f /D = 1.5, and f /D = 2. In order to demonstrate the effect of transverse 
expansion and laser divergence we show in Fig. 3a two sets of curves: one corresponding to f/D=l, the other to 
f/D=1.5. Each set consists of three curves: two curves with no target expansion and laser divergence, with /3g = 1 
and Pg < 1, and the third one with transverse expansion and Pg < 1. As one can see this effect significantly reduces 
the maximum ion energy by switching off the acceleration early. This switching off can clearly be seen in Fig. 3b, 
where the evolution of the reflected EM field amplitude at the foil is given for each of the curves in Fig. la. If the 
transverse expansion is taken into account, the EM field amplitude is much lower, which has a direct consequence in 
lower maximum ion energy. 

The utilization of an external guiding structure may relax the limits on maximum attainable ion energy. To model 
such interaction we solve Eq. assuming that the laser pulse is guided by a self-generated channel with a transverse 
size of wq = 0.9A {f/D = 1), wq = 1.35 {f/D = 1.5), and wq = 1.8 {f/D = 2), and the foil stays opaque for the 
pulse (|p| = 1). The solutions for such configuration are shown by dashed curves in Fig. 3c. One can see that the 
external guiding structure significantly enhances the maximum attainable ion energy, which is now limited by the laser 
group velocity in such a structure. This limit is illustrated in Fig. 3d, where the evolution of the reflected EM field 
amplitude is shown. The fields corresponding to the cases with guiding demonstrate the phase locking with the ions, 
i.e., they co-propagate with the foil, but are not able to transfer any energy to it. In principle, a composite target, 
consisting of a thin foil, followed by a near critical density (NCD) slab, may provide an example of such guiding. The 
laser pulse will accelerate the irradiated part of the foil in the self-generated channel in the NCD plasma [^. Though 
the foil density will drop due to the transverse expansion, the NCD plasma electrons being snowplowed by the pulse 
would provide an opaque density spike, which being pushed by the radiation pressure would drag the ions of the foil 
with it. Thus such configuration is similar to the one considered above: the laser pulse is guided with no diffraction 
and, although the density of ion decreases, the reflection coefficient is unity. 
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FIG. 4: The simulation setup: laser is focused on a thin foil, which is followed by a near critical density plasma. 


IV. COMPOSITE TARGET RPA: RESULTS OF 2D PIC SIMULATIONS 


In this section we present the results of 2D PIC simulations (using the code REMP (4l|), of an intense laser 
interacting with two types of targets: (i) a single thin solid density foil and (ii) a thin solid density foil followed by 
an NCD plasma slab (a composite target). As was shown in the previous section a composite target should provide 
guiding for the laser pulse and compensate for the transverse expansion of the foil, leading to higher proton energies 
than in the case of a single foil. The simulation box is 100 x 60A^, dx = dy = 0.025A, dt = 0.0125 x 2Trluj, and 
the number of particles per cell is 100. The Gaussian laser pulse is initialized at the left border with dimensionless 
potential qq = 100, waist w = AX and duration r = 30 fs, which corresponds to the average power of 1.8 PW. The 
pulse is focused with f-number f/D = 2at the left front size of the target, which is placed 16A away from the left 
border. The composite target consists of a fully ionized hydrogen foil and a hydrogen NCD plasma slab placed right 
behind the foil. The typical simulation setup is shown in Fig. [D 

In Fig. [5] we present the ion density distributions for the cases of a composite target and a single foil at different 
moments in time. The laser peak power is 1.8 PW and is focused with f-number=2 at the foil. The foil is 0.25A 
thick and the NCD plasma slab is 50A thick and its density is equal to critical one, Ue = Ucr- In both cases the 
transverse expansion of the target plays an important role, leading to a significant deformation of the foil and density 
redistribution. As a result the density near the laser axis is reduced. In the case of a single foil it means that the 
foil became transparent for radiation and the acceleration has stopped. In the case of a composite target it is not 
so, since the background electrons from the NCD plasma being snowplowed by the laser compensate for the density 
decrease. Such electron behavior is shown in Fig. [51 In the case of a single foil the electrons are subjected to transverse 
expansion, which is similar to that shown in Fig. [5] for the ion density distribution. In the case of a composite target, 
the electrons that were originally in the foil are also pushed to the side but the NCD electrons being snowplowed by 
the laser pulse form a dense shell in front of the laser, providing a robust relativistic mirror for proton acceleration 
up to the end of the NCD plasma slab. The different behavior of the electron component of plasma in the cases of a 
single foil and a composite target should translate into different values of the maximum ion energy. 

The evolution of the maximum ion energy in both cases is shown in Fig. [7]for two densities of the foil: rig = 225ncr 
(Fig. [7}i) and Ue = AOOricr (Fig. [Tb). One can see that in the case of a composite target the ion maximum energy 
is significantly higher, which is due to the longer acceleration length in the NCD plasma. In the case of a single foil, 
the target becomes transparent at f ~ 35, when, according to Fig. 5, the density on-axis drops several times below 
the initial one. The final energy for the single foil is in good agreement with the analytical results of the previous 
section, shown in Fig. 3. The maximum ion energy in the case of a composite target should be determined by the 
group velocity of the laser in the self-generated plasma channel. However, due to the fast depletion of the pulse as 
well as its reflection at laser-plasma interface, the laser group velocity can not unambiguously be determined from 
the results of PIC simulations. In this case, we chose the velocity of the of the laser-plasma interface, /3/, as the 
characteristic quantity for ion acceleration in the channel. Using results of Ref. we can find the normalized 
velocity laser-plasma interface from the condition that it takes the depletion energy time, tdepi, for the laser pulse tail 
to reach the laser-plasma interface. This yields 


i.e., 


Pi 




(32) 


li ~ 


Ig 

31/2 


(33) 


(for the definition of the group velocity of the laser pulse inside the self-generated channel in the NCD plasma, 7 ^, 
see Appendix A). For the case considered 7 / ~ 5, which is in good agreement with the results of PIC simulations. 
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FIG. 5: The evolution of the density of the ions originating from the foil during the laiser pulse interaction with a composite 
target (a-c) and a single foil (d-f) at t=35 (a,d), t=50 (b,e), t=75 (c,f). 
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FIG. 6: The evolution of the density of the electrons during the laser pulse interaction with a composite target (a-c) and a 
single foil (d-f) at t=35 (a,d), t=50 (b,e), t=75 (c,f). 


taking into account pulse reflection during the initial interaction with the foil. In Fig. [71 the energy 7 / — 1 (from PIC 
results) is shown by thin red curves. Note that 7 / — 1 remains constant during the pulse propagation through the NCD 
plasma, but after the pulse depletion becomes significant, it decreases, marking the end of the ion acceleration for 
both values of the foil density studied in simulations. We note that the ion energy is not able to reach the maximum 
value of 7 / — 1 due to the laser pulse depletion in NCD plasma. 

In order to further illustrate the differences between these two cases, we show in Fig. [ 8 ] the spectra of protons in 
the case of a composite target (Fig. |5^) and a single foil (Fig. [BJd) at different moments in time. While at t = 35 the 
spectra look similar, at t = 50 the difference becomes pronounced, since a single foil became transparent for radiation, 
marking the termination of the acceleration process. Moreover, at t = 75, in the case of a composite target, the 
protons gained almost 50% of their energy at / = 50, and, in the case of a single foil, the gain is as low as 10%. 

In Fig. we show that the protons can originate from different parts of a composite target. If we take such target 
apart, each of its parts, when irradiated by an intense laser, can generate high energy protons. We compared the 
maximum energy evolutions in all these cases to find that the protons from the foil, in the case of a composite target, 
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FIG. 7: The time evolution of the ion energy for a composite target (blue curves, 1) RPA and single foil (magenta curves, 2) 
RPA. The foil thickness is 0.25A with densities Ue = 225ncr- (a) and rie = lOOricr (b). The thickness of the NCD plasma slab 
is 50A and density is equal to ricr- The evolution of 7 / — 1 is shown by thin red curves (3). 




FIG. 8 : The evolution of the ion spectrum originating from the foil during the laser pulse interaction with a composite target 
(a) and a single foil (b) at t=35 (red curve), t=50 (blue curve), and t=75 (black curve). The parameters of the interaction are 
the same as in Fig. 3. 


acquire maximum energy of 2.4 GeV. Note that the energy of protons from a single NCD slab is in good agreement 
with the results of the analytical estimate (see Appendix B). Thus, the addition of an NCD plasma slab at the back 
of a thin foil, is a crucial feature of the target design. 

In Fig. [ni 3 we show the dependence of the maximum ion energy on the density of the NCD plasma. Since the 
acceleration is characterized by the velocity of the laser-plasma interface, it is plausible to assume that the energy 
of ions would scale as 77 — 1 scales with density. One can see from Fig. 8 b that maximum ion energy has similar 
dependence on density as 77 — 1. Here we took into account the fact that a significant part (~ 50%) of the laser pulse 
was reflected at the foil before the remaining part of the laser pushed the irradiated part of the foil inside the NCD 
plasma. However, for low densities, the similarity between 77 — f and the maximum proton energy fails, due to the fact 
that the laser propagation through such plasma can no longer be described as a propagation inside a self-generated 

channel [I3- 

Thus we identified the dependence of the maximum ion energy on the density of the NCD plasma and the fact that 
this energy is characterized by the velocity of the laser-plasma interface. We compared the cases of a single foil and 
a composite target, showing that in both cases the transverse expansion signihcantly modihes the target. In the case 
of a single foil, this modification leads to the termination of the acceleration due to the increased transparency of the 
target. In the case of a composite target the NCD electrons being snowplowed by the laser pulse form a dense shell 
in front of the laser, providing a robust relativistic mirror for proton acceleration. 















12 



FIG. 9: a) The evolution of the maximum ion energy in the case of (1) the protons from the NCD plasma slab (composite 
target), (2) the protons from the NCD plasma slab (no foil attached), (3) the protons from the foil (single foil target), and (4) 
the protons from the foil (composite target); b) The dependence of the maximum ion energy on the density of the NCD plasma 
slab for two values of the foil density ne = 225ncr (blue curve , 1) and rie = dOOricr (red curve, 2). The dashed curve (3) marks 
the maximum proton energy in the case of a single foil. The dotted curve is the dependence of 7 / — 1 on density. 


V. CONCLUSIONS 

We considered laser driven acceleration of ions in the radiation pressure acceleration (RPA) regime, taking into account 
the fact that the laser group velocity can be smaller than the light vacuum speed. We also included the deformation of 
the thin foil target, i.e., the transverse expansion in the equations of motion of the foil. This allowed us to identify two 
factors that limit the maximum attainable in the RPA ion energy: (i) the fundamental effect of the laser pulse group 
velocity being less than vacuum light speed, and (ii) the transverse expansion of the target, which plays a major role 
when tightly focused pulses are used. The first one results in the energy transfer from the pulse to the foil vanishing 
as the velocity of the foil approaches the group velocity of the laser. Thus the ion velocity can not exceed the laser 
group velocity. The transverse expansion of the target effectively decreases the density of the foil, eventually making 
it transparent for radiation. For realistic parameters of the laser foil interaction the effect of the transverse expansion 
dominates over the group velocity effect. 

We showed that the utilization of an external guiding may relax the constraints on maximum attainable ion 
energy. Such guiding should accommodate the transverse expansion of the ion component without the foil becoming 
transparent. We showed that a composite target, a thin foil followed by an NCD slab, is an example of such guiding 
structure. The NCD slab provided guiding of the laser pulse during the acceleration process as well as the electrons 
to replace the ones that are expelled due to the transverse expansion. The comparison of a single foil RPA and a 
composite target RPA shows that in the latter case the ions can gain the energy several times larger than in the 
former case, thus greatly increasing the effectiveness of the RPA regime of laser driven ion acceleration. In such a 
configuration the group velocity effects begin to dominate and determine the maximum achievable ion energy. 
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Appendix A: Laser propagation in Near Critical Density Plasma 

Here we address the propagation of an intense laser pulse in a near critical density plasma. The proposed composite 
target concept depends on the self-channeling of an intense laser pulse in such a plasma. The determination of the 
channel parameters, radius and length, the parameters of the laser pulse mode, which will propagate in this channel, 
are prerequisite for the study of the proposed mechanism. In principle, the propagation of an intense laser pulse in 
the relativistically underdense plasma can be approximated by the propagation of an EM wave in a waveguide [ 13 . 
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It is connected with the fact that, as the laser pulse propagates in such a plasma, the ponderomotive force of the 
laser expels electrons and ions from the region of the laser pulse. However due to much smaller charge to mass ratio 
the time of ion response is much larger than that of electrons, which leads to the formation of a positively charged 
region just after the front of the laser pulse. The expelled electrons are subject to the laser ponderomotive force and 
the attracting charge separation fields. This results in the formation of an electron density channel with high density 
walls. After the laser pulse passes, the ions also expand in the transverse direction forming an ion density channel. 
This is partly due to the response to the ponderomotive force and partly due to the coulomb explosion of a positively 
charged channel background. The laser radiation at the early stages of channel evolution is not able to penetrate 
these walls and is totally contained inside the channel. Then, the laser mode with lowest transverse frequency, an 
H-wave {E^ = 0), will propagate through this channel: 

= AJi{nr) cos {ujt-kx), = ^ 

where Ji is the Bessel function of the first kind and H = 1.84/i?; R is the radius of the channel, and A = 
(2H/a;) (mwa/e), here a is defined as the amplitude of the transverse electric field. The radius of the channel can be 
estimated from balancing the electron energy gains from the laser pulse EM field and from the field generated by 
the positive charge of ions in the channel: Rch = \/a {uj/ujp) X/tt. The group velocity of such an EM wave in the 
waveguide is /3g = [l — (1.84/2)^(wp/w)^a“^]. The total energy of such a pulse in the channel is 

= nR^Ta^mf.ncrK, (A2) 

where r is the pulse duration and the factor K comes from the integration over transverse coordinates and the duration 
of the laser pulse 


- Jo{^R)J2{^R)\ ^ (A3) 

Here Jo s-nd J 2 are Bessel functions of the zero and second order respectively. The relation between the amplitude of 
the laser pulse vector-potential and the laser pulse power, P = ££^/t, is 


_ / 2 P 


(A4) 


where Pc = 2m\le^ = 17 GW. For P = 1 PW and Ue/ncr = 1, the relation yields a ~ 116. As was shown in Ref. [l3| 
the depletion length of the laser pulse in the NCD plasma is Lch = a{ncr/ne)LpK, where Lp is the length of the laser 
pulse. Using Rch = \/a (uj/tOp) A/tt we can write the depletion length, normalized to pulse length, and the radius of 
the channel, normalized to laser wavelength, in terms of the laser pulse power: 


Ech 



Rch 


2l/3j^2l3 



1 fricrV^^ f 2P 

TT \ne ) \KPc 


K’ 

1/6 

(A 6 ) 


The 7 -factor, corresponding to the group velocity of the EM wave in the waveguide, neglecting pulse depletion and 
evolution, is 


7 


ch 

9 



(A7) 


If a thin foil is accelerated by the radiation pressure of an EM pulse inside this channel, then 7 °^ is the maximum 
achievable energy of the ions. 


Appendix B: Maximum ion energy via MVA mechanism 

In what follows we estimate the maximum ion energy that can be gained via the MVA mechanism in an NCD 
plasma [i^. In the MVA regime, the acceleration is achieved by an electric field at the back of the target, which is 
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generated by, and is of the order of, the magnetic field (Bch) of electrons, accelerated by the laser pulse in plasma in 
the forward direction, as they exit the target. This EM field, associated with the electron current leaving the target 
expands along the back surface of the target, while the magnetic field flux is conserved. The maximum energy gain 
of a charged particle in such an expanding field is She ~ ^ZHeBchRch- 

Let us assume that the electron current consists of all the electrons that initially were in the volume which is 

the volume of a cavity left in the wake of a laser pulse, ensuring that the fields of electrons and ions are compensated 
outside the cavity and the strong fields exist only inside. As the electrons move through the charged ion background, 
they experience a plasma lensing effect [d^, i.e., they are pinched towards the central axis and the radius of the electron 
beam is determined from the balance of the transverse electric field of the ion column, Ei = 27renei?b, which in the 
reference frame of an e-beam is E[ = ^eEi, and the self field of an electron beam, which is equal to A' = 2Tren'^Rt. 
From the condition E' = A' we obtain = Rch/le- Here 7 e is the Lorentz factor of the bulk of electrons accelerated 
forward. In the regime of laser pulse interaction with a NCD plasma, the electrons are continuously injected in the 
cavity behind the laser pulse, leading to the formation of an electron current, which is dominated by low energy 
electrons. 

The characteristic energy of these electrons is the injection energy, which can be obtained from the condition on the 
electron velocity to be equal to the group velocity of an EM pulse propagating in a waveguide of radius Ret'- 7 e = 7^^ 
(IAtI) . The magnetic field, generated by these electrons, is Bch{Rb) = ^ireneRchle at r = i?t,. Then the energy gain 
of an ion in the electric field at the back of the target is 

£He=m,2n^ZHe(^^^ (Bl) 

which gives an upper estimate for the ion energy gain in the MVA regime. For example, for a 1 PW laser pulse and 
Ue = Ucr, the protons energy is £h ~ 500 MeV. From Eq. (EH), we see that the ion energy scales with the laser pulse 
power as The scaling was derived assuming the optimal match between the laser pulse and the plasma, i.e., the 

target length is equal to the laser pulse depletion length. We should mention here that the MVA regime is relevant 
for ion acceleration from targets with densities not far from critical 0. For lower densities, laser pulse filamentation 
and Langmuir wave generation prevent efficient channel generation. At higher densities other mechanisms of laser ion 
acceleration come into play, and the target thickness, which is equal to the depletion length, is of the order the pulse 
length or smaller, which is outside the regime of MVA applicability. 
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